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, Abstract 
^ . 

^f) ' We consider pair production and decay of fundamental unstable particles 

in the framework of a modified perturbation theory treating resonant contribu- 
tions in the sense of distributions. Outcomes of a numerical simulation up to 
I the NNLO with taking into account universal massless-particles contributions 

are presented. 

^: 

r-| I In realistic descriptions of the processes of production and decays of fundamental 

unstable particles the gauge cancellations and high accuracy of the calculation of 
resonant contributions must be provided. A modified perturbation theory (MPT) 
[U 12] is one of the approaches intended for the solution of these problems. Its main 
feature is the direct expansion of the probability instead of the amplitude in powers 
^ I in the coupling constant, on the basis of distribution-theory methods. As the object 
to be expanded (the cross-section) is gauge invariant, the gauge cancellations in the 
^ ■ MPT must be automatically maintained. However the accuracy of the calculation 
O . of resonant contributions remains vague. For clearing up this question numerical 
simulations are required. 

In this report I present results of calculations in the MPT up to the NNLO in a 
^ ' model defined as the improved Born approximation for the process e~^e~ — )■ 7, Z — )■ 
■ — )■ W~^b W~b. The total cross-section of the process is 



a{s) = / -0 -;s // S1S2 a(s';si,S2) (l + ^c) , 



> 
O 




where 0(s'/ s; s) is the flux function, ct{s ; Si, S2) is the exclusive cross-section, 6c stands 
for soft massless-particles contributions, si and S2 are virtualities of the top quarks. 
With massless b quarks simin = S2min = M^r, Smin = 4M^. The cj) we consider 
in the leading log approximation. By 6c the Coulomb factor is meant, taken into 
consideration in the one-gluon approximation with conventional resummation, not 
affecting the Breit-Wigner (BW) factors (see details in |2]). 
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Table 1: The total cross-section in pb and in % with respect to the exact result. 
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0.6344 
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94.3% 
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0.2255 


0.1821 


0.2124 


0.2240(2) 




100% 


80.8% 


94.2% 


99.3(1)% 



In a{s ; Si, S2) the kinematic and BW factors can be extracted, 

a{s;si, S2) = ^6'(ys -^/ii -yi2)\A(^^^i^^ "^(-s; si, S2) p{si) p{s2) . (2) 

Further, $(s; Si, S2) niay be determined in the conventional perturbation theory, but 
the BW factors p(si)p(s2) must be determined as asymptotic series in the distributions 
sense in powers in the coupling constant a. For isolated p{s) in the case of a smooth 
weight, the series up to the NNLO looks as follows [Ij: 

2 

p(s) = 5{s-M^) + Taylor {p(s)}+ ^ c„ ^("^(s-M^) + 0(a=^) . (3) 

n = 

Here M is the top quark mass, 5(- ■ ■) is the 5-function, 5*^"^ is its n-th derivative. The 
^PV Taylor' means the naive Taylor expansion with principal- value prescription for 
the poles. Coefficients c„ within the NNLO are determined by three-loop self-energy 
contributions (for explicit expressions see [2]). 

Unfortunately, the weight in (|2]) is not smooth because of the 6'-function and 
the square root in the kinematic factor. However it can be analytically regularized 
by substitution [\{s, Si, S2)Y^'^ — > [\{s, Si, 82)]" ■ The singular integrals then can be 
analytically calculated and after removing the regularization the results become finite. 
Moreover, the expansion remains asymptotic 

The further numerical calculations, I carry out by means of rather general FOR- 
TRAN code with double precision. Parameters of the model I determine as follows: 
M = 175 GeV, Myy = 80.4 GeV. The results for the total cross-section above the 
threshold are shown in the figure (notation is obvious). The results in numerical form 
at typical ILC energies are presented in the table. 

It should be emphasized a very stable behavior of the NLO and NNLO approxi- 
mations at the energies above and near the maximum of the cross-section (this result 
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to a large extent is model- independent |3]). Near the maximum the discrepancy of 
the NNLO is within ±0.5%, which seemingly is sufficient for the ILC in the case of the 
top quarks. The further increase of the precision is possible at the proceeding to the 
NNNLO [2]. Another way for boosting the precision may be based on the combined 
use of the MPT, when the loop corrections are handled in the MPT framework but 
Born contribution in the conventional approach with the Dyson-resummation in the 
unstable-particles propagators (on the analogy of the actual practice of application 
of DPA). 
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